Overscreened Kondo fixed point in S = 1 spin liquid 
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We propose a possible realization of the overscreened Kondo impurity problem by a magnetic 
s = 1/2 impurity embedded in a two-dimensional S = 1 (7(1) spin liquid with a Fermi surface. 
This problem contains an interesting interplay between non-Fermi-liquid behavior induced by a 
U(l) gauge field coupled to fermions and a non- Fermi-liquid fixed point in the overscreened Kondo 
problem. Using a large-iV expansion together with an expansion in the dynamical exponent of 
the gauge field, we find that the coupling to the gauge field leads to weak but observable changes 
in the physical properties of the system at the overscreened Kondo fixed point. We discuss the 
extrapolation of this result to a physical case and argue that the realization of overscreened Kondo 
physics could lead to observations of effects due to gauge fields. 

PACS numbers: 71.27.+a, 71.10.Hf, 75.10.Kt, 72.10.Fk 



I. INTRODUCTION 

Impurity models constitute an important chapter in 
modern condensed matter physics. Since the original 
paper by Kondoi considering electron sea screening a 
single impurity spin, this problem has attracted signif- 
icant theoretical and experimental attention^— More 
recently, impurity physics has been studied in the con- 
text of strongly interacting systems. Numerous examples 
include^ 8 - an impurity in systems with vanishing density 
of state d 19 ' 20 high temperature superconductors^ and 
quantum magnets^ 2 - - — Quantum magnets are particu- 
larly versatile as a host system, having a large number 
of possible ground states with different low energy exci- 
tations. 

In this paper we consider a spin-half impurity em- 
bedded in a spin-1 quantum paramagnet with a spin 
liquid ground stated We consider the situation where 
the low energy excitations of the paramagnet are de- 
scribed by emergent fcrmionic excitations with a Fermi 
surface, coupled to a U(l) gauge field. This study is mo- 
tivated by the recent appearance of several 5=1 spin 
liquid candidate materials! 29 : 30 Theoretically, a number 
of spin liquid ground states for spin-1 system have been 
proposed^Ir— One possible scenario involves emergence 
of three fermionic excitations carrying spin-1 quantum 
numbers . 34 : 36 ' 38 Assuming that Fermi surfaces of these 
excitations are not destroyed by a pairing instability, we 
obtain the host system that is considered below. 

Impurity physics in a spin-1/2 spin liquids has been 
considered in the context of bosonic spin liquids^ 3 - al- 
gebraic spin liquids ; 24 ' 25 and spin liquids with a Fermi 
surfaced In particular, Ribeiro and one of us in Ref. l27l 
concluded that physics of a spin-1/2 impurity embedded 
in a spin liquid with 5 = 1/2 fermionic excitations is sim- 
ilar to that of the conventional Kondo problem^ In what 
follows we argue that a spin-1/2 impurity in a 5 = 1 spin 
liquid with a Fermi surface realizes overscreened Kondo 
physics. Although our results are qualitatively similar to 
the overscreened Kondo effect in conventional systems, 
there are observable differences due to the presence of an 



emergent gauge field coupled to spinons. 

Our findings suggest that an impurity in a 5 = 1 
spin liquid can be used to probe fermionic excitations. 
As these excitations do not carry charge, their exper- 
imental detection is a difficult problem. Different ex- 
perimental probes have been suggested in the context of 
spin-1/2 spin liquids i22r— We suggest that the realiza- 
tion of overscreened Kondo physics is a possible way to 
unravel physics of spin-one spin liquid, allowing probes 
of fermionic excitations, as well as the presence of an 
emergent gauge field. 

Overscreened Kondo physics is realized in multichannel 
Kondo models, where a single spin is coupled to N copies 
(flavors) of itinerant electrons^ On the one hand, such 
a generalization of original Kondo model may be seen as 
merely a theoretical tool, allowing a perturbative expan- 
sion in 1/N. On the other hand, the physics changes 
drastically depending on the interrelation between impu- 
rity spin length, s, and the number of flavors coupled to 
the impurity. When the number of flavors, N, is just 
enough or less than needed to screen the impurity spin, 
N < 2s, antifcrromagnetic coupling between the impu- 
rity and electrons flows to infinity in the infrared, mean- 
ing that at low temperatures impurity spin is screened to 
the maximum possible extent by electrons, resulting in 
Fermi-liquid behaviori 5 -^ However, in the overscreened 
regime, N > 2s, i.e. when there are more channels than 
required to screen the impurity spin, the system has a 
non- Fermi-liquid fixed point ,SiL2iII This state is charac- 
terized by singularities in different physical observables, 
such as impurity spin susceptibility, specific heat, etc. 
It is particularly interesting as a solvable example of 
a system with a non-Fermi-liquid fixed points Despite 
the rich and interesting physics, the overscreened regime 
of Kondo model has only few realizations (in particular 
quantum dots and two levels systems! 13 : 15 ' 16 ). Hence our 
system is also interesting as a possible implementation of 
overscreened Kondo physics. 

Qualitatively, the problem of a spin-half impurity 
hosted by isotropic 5 = 1 spin liquid looks similar to the 
conventional overscreened Kondo impurity model. When 
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coupled antiferromagneticaly, itinerant excitations car- 
rying spin-1 quantum numbers cannot screen the impu- 
rity. However, the presence of a gauge field effectively 
enforcing a single occupancy constraint for fermionic ex- 
citations, makes these two problems different. Even 
without the impurity, fermions are in a non-Fermi-liquid 
regime^— due to the gauge field. The fermion prop- 
agator is dressed by a singular self-energy, so there are 
no well defined quasiparticle excitations in the system. 
This is manifested, for example, in the singular behavior 
of the specific heat C cx T 2 / 3 in two dimensions at low 
temperatures , 45 1 49 

Coupling the impurity to fermions with non-Fermi- 
liquid behavior allows us to study the interplay between 
the gauge field induced non-Fermi-liquid behavior and 
the Kondo non-Fermi-liquid fixed point. The conven- 
tional approach to the Kondo problem is either an exact 
solution by mapping it onto one-dimensional problcm^iii 
or 1/N expansion. Both methods are not directly appli- 
cable in our case. The presence of gauge field impedes 
the mapping of our model to a one dimensional problem 
in the radial channel. On the other hand, a rigorous 1/N 
expansion is not possible, due to singular self-energy cor- 
rection a 50 ' 52 . The latter issue has been recently resolved 
in the paper by Mross et. al.^ where a controlled dou- 
ble expansion scheme has been provided. It combines 
the 1/N expansion with an expansion in another small 
parameter (related to the dynamical critical exponent of 
the gauge field). 

We adopt the recently developed double expansion 
method^ to our problem. Since the double expansion 
includes the large N limit, we expect to have a perturba- 
tively accessible fixed point. At leading order, the gauge 
field does not affect the position of this non-Fermi-liquid 
Kondo fixed point. However, it leads to corrections to 
the scaling dimension of the Kondo coupling. Assuming 
that the results obtained using the double expansion in- 
terpolate to the physical case, we conclude that physical 
properties such as impurity spin susceptibility, specific 
heat, etc. are still characterized by singular behavior. 
Unlike the case of the Kondo model in the regime of per- 
fect screening j2I where the coupling to the gauge field 
has no consequences to leading order in 1/N, in our case 
the gauge field influences Kondo physics. 

The rest of the paper is organized as follows. In the 
remainder of this Section, we introduce the basics of 
our model, diagram technique and briefly explain the 
idea behind double expansion. In Section [H] we first re- 
view known results for the /3-function in the overscreened 
Kondo problem without the gauge field. Afterwards, we 
calculate the /3-function with the gauge field and study 
the changes in scaling behavior of different physical quan- 
tities. Finally, in Scction lHll we discuss the extrapolation 
of our findings beyond the double expansion, and com- 
ment on possible experimental realizations and experi- 
ments to detect Kondo physics. Details regarding the 
calculation of corrections to the /3-function due to gauge 
field are given in the Appendix. 



A. Spin liquid with fermionic excitations and 
impurity 

Our starting point is a spin Hamiltonian on a lattice 
consisting of spin-1 sites, 

# sp in = ]T [Jij Si ■ Sj + Ka (Si ■ S,-) 2 ] + • • • , (1) 

where ellipses denote other possible terms needed for sta- 
bilizing the U(l) spin liquid phase with a Fermi surface. 
We do not specify the lattice, since only an effective low 
energy description of a spin liquid phase is relevant in 
what follows. However we note that recent work 3 -^ shows 
evidence for the possibility of stabilizing such a phase on 
a triangular lattice with nearest neighbors bilinear, and 
biquadratic spin interactions along with ring exchange 
terms. In the large-t/ limit of a half-filled Hubbard model 
one can easily derive the coupling between a given lattice 
spin at site i and an impurity spin in a form 

Hi ni p = JfcSi ■ S; m p. (2) 

This procedure leads to Jk > corresponding to anti- 
ferromagnetic coupling. 

We assume that Hamiltonian in Eq. (fTJ) has a spin liq- 
uid ground state with fermionic excitations and a Fermi 
surface. The low energy description of such a state is a 
theory of fermions strongly coupled to a U(l) gauge field. 
Referring the reader to the literature for more detailed 
discussions^ we only summarize results. A spin-1 op- 
erator at a _given site is represented using three fermion 
operators, fx, A = 1, 2, 3 as 

3 

*= E fk iXp f*p> ( 3 ) 

x,p=l 

with I Xp being the set of three spin-1 matrices (gener- 
ators of SU(2) in spin-1 representation)^ In order to 
remove unphysical states from the Hilbert space, intro- 
duced by the representation in Eq. §5§ , one has to enforce 
a single occupancy constraint on each site. Fermionic fx 
are the low energy excitations of the spin liquid, carrying 
spin-1 quantum numbers. In addition, the low energy de- 
scription contains a U(l) gauge field, coupled to fermions 
fix and enforces the single occupancy constraint. 

Before proceeding further, let us reiterate the question 
of interest. We want to understand if the non-Fcrmi- 
liquid fixed point of a conventional overscreened Kondo 
model is changed by the presence of the gauge field in 
the bulk. The model outlined above provides us with a 
particular setup to study the influence of the non- Fermi- 
liquid bulk on the overscreened Kondo fixed point. How- 
ever, in order to have control over calculations we need to 
resort to the large- A limit. The crucial requirement for 
the generalization procedure is to retain the presence of 
the overscreened Kondo fixed point. We choose a model 
with N species of spin-half fermions, fi am with a = t, 4-, 



and m = 1 ... N as a large- N generalization. This is the 
simplest model which allows for controllable calculations. 

The corresponding Lagrangian for our generalized 
model may be split into a fermionic part (including cou- 
pling to gauge field and impurity spin) , and a gauge field 
Lagrangian, 

L -^fcrmion ~t~ -^gaugc- (4) 

The generalized fcrmion Lagrangian becomes: 



L 



fcrmion 



/ dr E (/* 

(d T -e k )fk C 



h . m . a 



E^+fc™^) ■ «<?/fe-fam - -JfS(0) ■ S 



(5) 



where we use imaginary time. In accordance with the 
discussion above, fermion operators fi am now carry spin- 
1/2 quantum numbers (a = til)- We omit the time 
component of gauge field from the coupling, since it is 
screened^ and do not write the diamagnctic term, in- 
cluding it in the gauge field propagator. The fermion 
spin at r = is 



W) = jr E fL 



p/3m i 



(6) 



with er = (a x ,<j v ,(t z ) being the set of three Pauli ma- 
trices, and Af being a number of sites in the lattice. 
In what follows, Greek indices label spin projection, 
a, j3, . . . = t; 4-) whereas Latin indices m,n,... = I . . . N 
label channels. 

The gauge field Lagrangian is 



1 f dqdus 



a* (<?, ^D- 1 (q, iu))aj (q, u>), (7) 



where the time component of the gauge field is excluded. 
The bare gauge field propagator is zero, since the gauge 
field is not dynamical but rather represents fluctuations 
around the mean field ansatz. However, nontrivial dy- 
namics are generated if one accounts for coupling to 
fermions, leading to a non-zero D'^ 1 (q,iui), discussed in 
Section DZ3 



B. Diagram technique 

The impurity spin is conveniently represented via 
fermionic operators, 



E c «^r c /3' 



(8) 



a, 13 



where cl ^ {cf,i) are creation (annihilation) operators of 
spin up or down pseudofermions <2 In what follows we use 



G(k, iuj) 



F(i 



8 7 



D(q,iu) */\/\/\/\/\. 



a 

< 



FIG. 1. Summary of rules for diagram technique. Solid, 
dashed and wavy lines represent fermion, pseudofermion and 
gauge field propagator respectively. Also, interaction vertices 
of fermions with gauge field, X a g, and fermions with impurity 
pseudofermions, r a a 7 a, are shown. All objects are diagonal 
in flavor indices, which are thus suppressed. 



the term "pseudofermions" to distinguish the operators 
c a from the operators fkam, which describe low energy 
excitations in the spin liquid. A faithful representation 
of spin via fcrmion operators requires an additional con- 
straint to exclude doubly occupied and empty states from 
the Hilbcrt space. However, in the case of a single spin- 
1/2 operator, s, in Eq. ©, gives zero when acting on 
unphysical states in the Hilbert space. Therefore, one 
can ignore the constraint in this case^ writing the impu- 
rity Lagrangian as 



-^imp — / C m (c^ T 



A^imp)c?n , 



(9) 



where fj,i mp = +0 is the chemical potential for impurity 
pseudofermions . 

The rules for diagram technique, following from 
Eqs. ©-([9]) are summarized in Fig. Q] Propagators for 
fermions and pseudofermions along with interaction ver- 
tices are given by 



$ a (3$ win 



F a p(iw) 



1 a/37<S 



iuj — — S(iw) ' 



iuj — ' 
Jk 



V h 5 a gS, 



{3 "mri j 



(10a) 

(10b) 

(10c) 
(lOd) 



where = — is the fermion energy relative to 
the Fermi surface. The self-energy, included in fermion 
Greens function [Eq. (|10a[) ] is discussed below. We note 
that interaction between fermions and the impurity is 
local in real space. Therefore, in Fourier space, the mo- 
menta of two fermion operators in the impurity interac- 
tion vertex [Eq. (|10c[) ] are unrelated. Fermion propaga- 
tor and interaction vertices are diagonal in flavor indices. 
Thus the only contribution of flavor indices is an extra 
factor N for every loop of fermions, and they will be 
suppressed in the remainder of the paper. 
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(a) 




SL 



function. This spoils naive power counting in the 1/N 
expansio n 50 ' 52 requiring a summation of an infinite se- 
ries of diagrams of a particular topology (genus) at any 
given order in 1/N. However, if we assume a gauge field 



dynamical exponent 
scaling limits 



54 



lb 



2 + e, and take the double 




(b) 



FIG. 2. (a) Self-energy of the gauge field due to interaction 
with fermions. Second diagram describes diamagnetic contri- 
bution, (b) Self-energy of fermions due to interactions with 
gauge field in the leading order in 1/iV. 



0, 



TV ->• oo, 



sN = const, 



(13) 



we obtain finite A cx l/(Ns) in Eq. (jI2b[) . rather than 
A cx 1/N —> 0. The absence of the factor 1/N in front 
of the self-energy restores naive power counting, where 
the gauge field interaction vertex contributes I / \/N and 
each fermion loop gives a factor of N. 

Finally, before proceeding further, we rewrite H(iu}) in 
a simplified form, valid in the double scaling limit, 



C. Double expansion 

We briefly review the double expansion framework in- 
troduced in Ref . [HTl First, we specify dynamically gener- 
ated propagator of the gauge field. To leading order, the 
propagator is given by the fermion bubble with current 
vertices along with diamagnetic term shown in Fig. [5] (a). 
In the Coulomb gauge, Va = 0, the propagator is trans- 
verse and can be written a s 46 ' 48 ' 51 

D^iaM = Uu ^jf\ D^(q,uj), (11a) 

^ - 1 (g, lW )= 7 M +Xo ^-i j (ii b ) 

with zi, = 3 and 7 = 2n/kp, xo = l/(24-7rm) for fermions 
with quadratic dispersion. Note that the Landau damp- 
ing term is non-zero only for \u\ < VfQ- We assume Zb 
takes general value 2% < 3 and use it as a control param- 
eter. This approach is consistent because terms q Zb ~ l for 
Zb < 3 are non-local. Since Zb is not going to be renor- 
malized within perturbation theory, it is a valid control 
parameter. 

The singular form of the gauge propagator [Eq. (fTTT) ] 
leads to a singular self-energy correction for fermions. In 
the leading order in 1/N, the diagram in Fig. H](b) gives 



us 



48.50.51 



= -i\\uj\ 2/Zb signw, 



e 2 v F 



N 7 47rsin^ \xo 



2/zf, 



(12a) 
(12b) 



For Zb > 2, the self-energy is more important than 
the bare itu term in Greens function [Eq. (|10aj) ]. when 
|cj| < ujq. The energy scale luq is set by a combination 
parameters 7, xo and vf [see Eq. (fT5j)] and is of order of 
Fermi energy, the only energy scale related to fermions. 

When the self-energy, Eq. (fl~2"j) , is singular compared 
to the bare frequency dependence of the fermions' Greens 
function, a factor of 1 /N in the fermion self-energy leads 
to an extra power of N in the numerator of the Greens 



y ' Ne 



^0 



e/2 



where scale coq is explicitly given by 



/ 2 \ 2/e 
_X0( e V F X 

° 7 V 271 " 2 Xo, 



(14) 



(15) 



II. PERTURB ATIVELY ACCESSIBLE FIXED 
POINT 

The renormalization group (RG) approach in conjunc- 
tion with 1/N expansion has been proven to be fruitful 
when applied to the conventional Kondo impurity prob- 



lci 



,2-4.11.55 



The renormalization procedure is defined 



with respect to the fermion bandwidth, D. Eliminating 
states far away from the Fermi surface, one studies the 
induced flow in the dimcnsionless coupling g = vJk (y is 
the density of states assumed to be constant within the 
whole band). The corresponding fi- function is defined as 



d\ogg 
dlogD 1 



(16) 



and can be calculated perturbatively in g. This simplifi- 
cation is brought by the 1/N expansion and is justified 
in vicinity of fixed point located at small g* cx 0(1/N). 

When there is a gauge field coupled to fermions, the 
RG approach still can be applied. However, it requires 
some modifications. The usual 1/N expansion has to be 
replaced by the double expansion discussed above, but 
the RG flow is still defined with respect to bandwidth, 
D. Coupling of the fermions to the gauge field, e 2 , is 
treated as a constant, since a single impurity can not 
change its flow under RG. Likewise in the conventional 
Kondo problem, there exists a perturbatively accessible 
fixed point. After briefly reviewing the calculation for 
RG flow in conventional Kondo problem, we calculate the 
^-function in the presence of the gauge field and obtain 
physical properties in the vicinity of the fixed point. 
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FIG. 3. Diagrams contributing to the /3-function in the lead- 
ing order in 1/N. Diagrams (a) and (b) describe corrections 
to the vertex in the second and third orders of perturbation 
theory (symmetric counterpart of diagram (a) with direction 
of one of the fermion line changed is not shown) . Diagram (c) 
is the correction to the self-energy of pseudofermions, con- 
tributing to /3-function via renormalization of Z-factor. 




FIG. 4. Two types of vertex corrections in the leading order 
in 1/N due to gauge field, (a) Example of vanishing diagrams 
with a single gauge propagator connected by at least one end 
to the internal line, (b) Non-vanishing corrections, represent- 
ing a new non-local vertex (first diagram) and example of 
diagram leading to its renormalization (second diagram). 



A. /3-function in conventional Kondo problem 

While reviewing the RG procedure for the usual Kondo 
impurity problem we mostly follow Refs. [U and[^. Renor- 
malization of the dimensionlcss coupling g in the leading 
order is given by diagrams shown in Fig. [3J Diagrams 

(a) and (b) in Fig. [3J represent corrections to the bare 
interaction vertex in the second and third orders of per- 
turbation theory. These are the only diagrams up to the 
third order, which are logarithmically divergent and thus 
renormalize the coupling. Note, that diagrams (a) and 

(b) describe the contributions of the same order, since 
the latter diagram in addition to extra power of g oc 1 /N 
has a factor of N from the fermion loop. Diagram (c) 
in Fig. [3] describes renormalization of Z- factor of pseud- 
ofermions and also contributes to the /3-function. 

Calculation of the diagrams in Fig. |3j gives the (3- 
function: 



dloe 



N 



dlogD 



9 3 + 



(17) 



where ellipses denote terms 0(l/N 3 ) coming from higher 
order diagrams. These terms have the form C\Ng A + 
C2 A 2 g 5 with coefficients readily available in the lit- 
eratu re 6 ' 11 ^ 56 and listed in Appendix [5] One can eas- 
ily solve for a non-trivial unstable fixed point of this f3- 
function: 



9 =N + ---> 

Ao = -W) = | + 



(18) 
(19) 



where Ao is the negative slope of the (3 function at this 
fixed point. Ellipses here stand for terms 0(l/N 2 ). We 
see that g* is indeed small in 1/N, justifying the use of 
perturbation theory. 



B. Correction to /3-function due to gauge field 

As we shall see, within the double expansion frame- 
work, the gauge field produces a small correction to the 
regular /3-function. Therefore, it suffices to consider the 
effect of the gauge field to leading order. 

There are two types of effects related to the gauge field. 
First, the gauge field destroys well-defined quasiparticle, 
leading to non-Fermi-liquid behavior. This is manifested 
by the singular self-energy due to the gauge field in the 
fermion propagator [Eq. ([13] . Therefore, one has to re- 
calculate diagrams in Fig. [3J using the fermion propagator 
which contains the self-energy. A lengthy but straight- 
forward calculation (see Appendix [S] for details) yields 
an answer identical to the case without gauge field, how- 
ever, with a modified divergent logarithm. Namely, the 
standard log-divergent contributions are replaced by 



log — 



log 



D 



1 



2N 1 + {Ne)- 



(20) 



where energy scale ujq cx D was defined in Eq. (|15[) . 

Another effect of the gauge field is the appearance of 
vertex corrections. All diagrams describing vertex cor- 
rections can be split into two classes with representatives 
of each class depicted in Fig. [4] (a) and (b) respectively. 
Diagrams belonging to the first class have at least one of 
the ends of the gauge field propagator connected to the 
internal fermion line. In Appendix [B] we show that due 
to the transverse character of the gauge field propagator 
and the locality of interaction with the impurity, all dia- 
grams of this type with single gauge propagator exactly 
vanish^ 

In the vertex correction diagrams attributed to the 
second class, the gauge field propagator connects two 
external lines. One can think about these diagrams 
as describing a new interaction vertex [first diagram in 
Fig. U (b)] and its renormalization [all other diagrams of 
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this type]. This new vertex contains an additional small 
factor 1/N, compared to the original impurity interaction 
vertex. However, it is non-local, since it depends strongly 
on the relation between outgoing and incoming fermion 
momenta, k and p respectively. The vertex is logarithmi- 
cally divergent when the transferred momentum \k — p\ 
is close to 2fcp and is small otherwise! 46 ' 48 ' 51 The flow 
of this vertex to leading order in 1/N is identical to the 
flow of the standard vertex. Thus it does not influence 
the scaling in the vicinity of the fixed point. The effect of 
this vertex is to provide subleading corrections to differ- 
ent observables (due to the extra factor 1/N). Therefore, 
in what follows we do not consider this vertex. 

As we demonstrated, no new diagrams contribute to 
the /3-function up to order 1/N 3 . Calculation of dia- 
grams in Fig. [3] with self-energy included into fermionic 
propagator gives us 



00;) = (!-«) [g 



N 



(21) 



where ellipses represent terms 0(1/N 3 ). These terms 
have to be included for consistency, since k cx but 
are identical to those in the /3-function without gauge 
field, Eq. p7|) [account for gauge field in these terms will 
produce corrections 0(1/N 4 )}. 

The obtained /3-function, Eq. (|2"Tj). differs from the 
/3-function without gauge field, Eq. (fT7|) . by terms 
0(1/N 3 ). This correction does not shift the fixed point 
g* even at order 1/N 2 , compared to fixed point Eq. (|T8|) . 
However, it modifies the slope of /3(g) at the fixed point, 



Ai = — (1-k) + ... = (1-k)A , 



(22) 



compared to the slope for the case of the conventional 
Kondo problem, A [Eq. (|f 91) ]. Below, the slope of the 
/3-function, Ai will be used to determine the flow of 
the coupling in the vicinity of the fixed point, as well 
as the singular behavior of different measurable quanti- 
ties^ Therefore, the difference between Ao and A mod- 
ifies the behavior of different observables compared to 
conventional Kondo problem. 



C. Observables 

In order to understand how the non-Fermi-liquid fixed 
point manifests itself in observables, we first find the 
dependence of the running coupling constant gn(uj) 
on u. It can be determined from the flow equation 
dgn(uj)/dloguj = —(3 (go), by employing results for j3- 
f unction and its slope at the fixed point. Denoting the 
bare value of coupling at to = D, as gj^(D) = g, we 
have^i 



9r(u) = g 



A 



(23) 



where A is the slope of the /3-function which depends 
on the presence of the coupling to the gauge field. The 



position of the fixed point, g* , Eq. (fT8|) is not influenced 
by the gauge field. The Kondo temperature is TV = 
Dg N / 2 exp(-l/g) and C = (g*) 1 + NA / 2 exp(-A/g*). Wc 
assumed here that to < Tk and that the initial value of 
coupling, g, is small. 

The power law behavior of the running coupling leads 
to a similar behavior in different physical quantities. Cor- 
rections to different measurable properties within pertur- 
bation theory can be expressed as a series in coupling 
g. Applying the renormalization group to this series re- 
sults in singular behavior as a function of frequency or 
temperature with exponent proportional to A. For the 
case when there is no gauge field present, this program 
has been implemented in Refs.fTll. Generalization of this 
procedure to the case with the gauge field is straightfor- 
ward. 

The main effect of the gauge field is always related to 
the different values of slope of /3-function, A. Without 
a gauge field A = Ao is given by Eq. (flU)) . When there 
is a gauge field, we have A = Ai, specified in Eq. (|2"2"j) . 
While for thcrmodynamical quantities this is the only 
effect, transport properties and other quantities acquire 
small corrections to prefactors which are not given here. 

Calculating the contribution of impurity to the imag- 
inary part of self-energy of fermions gives the scattering 
rate due to the impurity. As a function of frequency, it 
acquires a cusp at oj = 0, j/t _1 (w) oc I — N((lo/Tk) a (cf. 
with a Lorentzian shape for a Fermi liquid fixed point). 
The correction to the resistivity due to Kondo interac- 
tion has a similar form, however, it is of little interest 
due to neutral character of fermionic excitations in spin 
liquid. The correction to the heat conductivity is poten- 
tially more interesting, 



5 k 



th 



A 



(24) 



as it can be distinguished using its proportionality to the 
impurity concentration, nj mp . 

Also one can calculate corrections to different ther- 
modynamic properties. A rigorous calculations of the 
self-energy allows us to find impurity specific heat with 
a critical exponent a = 2 A: 



a, 



2 S \T K 



2A 



(25) 



Magnetic properties, such as the impurity susceptibility 
as temperature T — > and dependence of magnetization 
on the field h = (j,bH at T = are given by 

2A 

- ) , (26) 



N( 



M = — - — 



h 

f K 



1 ( T 
T \Tk~ 

A 



(27) 



Likewise, it is possible to find an expression for fcrmion- 
fermion, fermion-impurity and impurity-impurity suscep- 
tibiliticsi 1 -. Lastly, we list results for Ximp( w ! T)/u> which 



7 



is a contribution to the NMR relaxation rate due to 
the impurity. Its behavior is again specified by A, and 
forw«T 



imp 



— oc T 2A - 2 . 



(28) 



III. DISCUSSION 



We have investigated the effect of the gauge field 
strongly coupled to fermions at a non- Fermi liquid over- 
screened Kondo fixed point. Using the double expan- 
sion framework, we demonstrated that the gauge field 
does not alter the position of the perturbatively acces- 
sible non- Fermi-liquid fixed point, but leads to correc- 
tions to exponents characterizing the behavior of different 
physical properties in the vicinity of the fixed point. In 
particular, it "softens" the non-analytic behavior of spe- 
cific heat, magnetization, spin susceptibility, compared 
to those for a Kondo problem without the gauge field. 
The physical origin of this effect is the "smearing" of the 
sharp quasiparticles by the gauge field. 

Let us discuss the extrapolation to the physical case. 
In order to have a control over calculations, we worked 
in the double expansion limit, Eq. (|13[) . with N species 
of spin-half fermions. We note, that if the coupling to 
gauge field was absent, the considered model for N = 
4 corresponds to the one channel of spin-one itinerant 
moments coupled to impurity.—^ The same equivalence 
was checked to hold in our perturbative calculations of 
/^-function for the case when there is a coupling to the 
gauge field. 

Thus, we expect that the physical case corresponds to 
N = 4, e = 1. Assuming that our results can be extrap- 
olated to these values of N and e, we can argue that the 
non-Fermi-liquid Kondo fixed point is not destroyed by 
the presence of a gauge field. However, we expect sin- 
gularities in different physical properties related to the 
non-Fermi-liquid fixed point to be weakened compared 
to their values without gauge field. In such a case, the 
realization of overscrccned Kondo physics in S = 1 spin 
liquid may be used not only to observe neutral fcrmionic 
excitations, but as evidence for the presence of a gauge 
field. Indeed, non-Fermi-liquid behavior may be used as 
an indication of fermionic excitations present in the sys- 
tem. At the same time, the difference of observed scalings 
from those for the case without a gauge field&i 2 - may be 
used as a litmus test for the presence of a gauge field cou- 
pled to fermions. From an experimental point of view, 
specific heat (proportional to impurity concentration) , as 
well as spin susceptibility and NRM relaxation rate are 
the most promising probes. 

It is instructive to compare the role of the gauge field in 
our case to the case of the Kondo model in the regime of 
perfect screening, Ref. H3- In the latter case, the system 
flows to the infinite coupling fixed point, and the results 
of Ref. show no changes in impurity specific heat and 
spin susceptibility due to the presence of the gauge field. 



Finally, we discuss possible experimental realizations 
of our proposal. In a recent experiment o 29 ' 30 materials 
that could possibly realize the spin liquid with fermionic 
excitation a 34 ' 36 ' 38 has been found. One can speculate on 
the possible stabilization of U(l) spin liquid phase in the 
same or similar type of materials. The presence of spin- 
half impurities in such a phase would realize the scenario 
considered in our work. Another way to implement the 
discussed physics is to go to lower dimensions. A gap- 
less phase for spin-1 with bilinear and biquadratic inter- 
action has been established for a certain range of cou- 
plingsj^r— . A spin-half impurity in such a chain is ex- 
pected to realize overscreened Kondo physics. A detailed 
consideration of this problem will be presented elsewhere. 

M.S. acknowledges discussions with David Mross and 
John McGreevy. This work was supported by NSF DMR 
1005434 (T. S.) and NSF DMR 1104498 (P. A. L.). 



Appendix A: Calculation of diagrams for /3-function 

In this Appendix we present calculation of diagrams 
in Fig. |3] with fcrmion propagator, Eq. (|10a[) , contain- 
ing self-energy due to gauge field. Detailed calculation 
of these diagrams without gauge field can be found, for 
example, in Refs. d and|||. 

First we consider diagram in Fig. [3] (a), describing 
second order correction to the dimcnsionless coupling 
g = vJk- We will be interested only in the logarithmi- 
cally divergent part of the diagram. Using zero tempera- 
ture Matsubara diagram technique and implying summa- 
tion over repeated indices we can write for the correction 
to impurity interaction vertex: 



-,(al) 
a/3715 



Jk 



(cr ap ■ (T 7A )(<T p ^ ■ (T\s) 

— — T F(-iu> 1 )G(k,iu> + iu> 1 ). (Al) 

(Z7T ) 



Symmetric counterpart of diagram (a) with flipped direc- 
tion of the propagation of the pseudofermions (not shown 
in Fig. [3]) gives us, 



^(a2) 
' a/3j5 



^^F(iwi)G(fc,m> + ^i). (A2) 



(27T) 



After integrating over ui\ and changing integration vari- 
able from k to £ = e*. — /x, we have similar expressions 
for both diagrams: 



-,(al,2) 



Jk 



(=F2<t Q( 3 • cr^s + 3<5 a/ 3<5 7l 5) 



x v I d£ 



*(±0 



w /2 )-e' 



(A3) 
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We perform integration over £, retaining only logarithmi- 
cal part. Collecting results for both diagrams and going 
to real frequency domain, we get: 



-Au 



Jk 



(<T a ■ CTjs) 



X loo 



D 



1-1(1 + w-,\^\ e ' 2 ) 



(A4) 



We expand logarithm in e to the leading order and collect 
both terms into single logarithm again: 



-hi ,/c2 



/dco 
-^G{k 1 ,iuj 1 )G(k 2 ,iuji - iw)0(-£ fe2 



(All) 

To simplify further calculations, we expand in e and 1/N. 
Self-energy, Eq. (fl"4)) . expanded to the leading order in e 
becomes: 



(A12) 



Inserting this into fermion Greens function, Eq. (|10a|) . 
and expanding in 1/iV we get: 



log 



D D oj 

1 \ui \ E / Z \ \u>\ M 



D 



where K is small, k oc 0(1/N), 

1 1 



2W 1 + W' 



= log — | (A5) 



(A6) 



Finally, we have 



D 



(A7) 



where bare vertex r Q( g 7< 5 is defined in Eq. (|10c|) . and we 
retained only logarithmically divergent terms. Alterna- 
tively, we could expand in Greens function in 1 /N already 
in Eqs. (|A1[) - (|A2[) . reproducing the same result. 

Calculations of vertex and Z-factor renormalization, 
described correspondingly by diagrams (b) and (c) in 
Fig. [3] are very similar. Indeed, in order to get impu- 
rity pseudofcrmions Z-factor, Z lmp , we have to differen- 
tiate self-energy over w, what may be thought of as an 
introduction of additional vertex with zero incoming fre- 
quency. Therefore, below we present only details on the 
calculation for the derivative of self-energy and list the 
result for the vertex renormalization. 

Correction to the impurity self-energy described by di- 
agram Fig. [3] (c) is written as 



E imp (iu;) = -6N 



Jk 



dki dwi dk 2 dQ 2 



A J J (2tt) 3 (2tt) 3 
x G(ki,iui)G(k 2 ,iuJi + iCl 2 )F(iuj + ifl 2 ), (A8) 

where we omitted spin indices of external pseudofcrmions 
and associated <5-function. Renormalization of Z lmp is 
given by the derivative of self-energy, 

d£ imp (» 



&Z lmp 

o{iu) 

Integrating over in Eq. (|A8[) . we have 



(A9) 



5Z-, 



imp 



Jk \ d f dh\ dk,2 



(A10) 



G(ki,iuji) = G{ki,iuj{) 
1 



2N 



iu>i log 



where G{ki 1 iuj{) is defined as: 
G(ki,iuji) 



^i(! + lfe)-Cfei 



[G(fc 1;lWl )] 2 , (A13) 



(A14) 



Finally, expansion of the product of Greens function in 
Ifei,fe 2 ,^> E q- 1 A11 I> - g ive s us: 

/duj 
-^—0(-^k 2 )G(k 1 ,iu} 1 )G(k2,ioJi - iui) 



1 iwi log 

2A 1 B 



-— i(u)t -w)log 



w 



UJ\ — UJ 



G(k ll iuj 1 ) 
G(k 2 , icoi — ilo) 



(A15) 



After integration over zeroth order term in (|A15j) 
yields 



r (0) 

~"k\ ,/C2 ,ioj 



(A16) 

This is inserted into Eq. (|A10[) . After integration over 
momenta extra factors (1 + -^-) drop out and we re- 
produce the answer for the case without gauge field, 
5zP> p = Wg*/8i 0g (D/\ U \). 

Frequency integration for terms proportional to 1 /N in 
Eq. (|A15[) results into cumbersome expression. However, 
after integrations over ^ and and extracting log- 
divergent part we obtain 5Z^J p = —3nNg 2 /8 log(cj /M), 
where k is defined in Eq. (|A6|) . Combining SZ^ p and 
SZ^ p , we have for impurity pseudofcrmions Z-factor: 



Z imp = 1 + SZ imp = 1 + -Ng 2 log 



D 



ii— «,,,« ' 



(A17) 



Correction to the impurity interaction vertex, diagram 
Fig. |3] (b) is calculated in a similar way. Resulting con- 
tribution to the interaction vertex is 



D r 



(A18) 
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v(k) 



k + q,U! + Q 



fc, 



FIG. 5. Part of diagram with vertex corrections that makes 
the diagram vanish. 



Finally, renormalized coupling is 

9 + Sg 



9R 



Z 



imp 



(A19) 



where Z lmp is given by Eq. (|A17|) . and Sg can be read 
from Eqs. and $MB- 



exact to the order 1/N 3 : corrections to subleading terms 
from to the gauge field are of order 0{\/N A ) and thus 
can be ignored. 



Appendix B: Vertex corrections 

In this Appendix we demonstrate that a subset of ver- 
tex corrections where gauge field propagator is connected 
to internal fcrmion Greens function vanish. Two exam- 
ples of such diagrams are shown in Fig. [4] (a). It suffices to 
consider a part present in all diagrams, consisting of two 
Greens functions and a single gauge field vertex, Fig. [5] 
Using notations adopted in Fig. [5j we can write for the 
integral over momentum k 



^=( 9 -^Alo gn ^- K . (A20) 

Using that loq oc D, we obtain the /3-function: 
rflogg / 2 N 3 

N N 2 «. 

~-(l + log2) g 4 + x3 5 , (A21) 

where in the second line we included subleading terms 
0(1/A 3 ) obtained in Ref. GH. Note, that Eq. ([Mil is 



dk x dk v v y (k)G(k, iuj)G(k + qe x , iuj + ifl), (Bl) 



where we assumed that q has only a;-component, q \\ e x 
and used fact that gauge field is transverse. Note that 
integration over k docs not involve any other functions 
due to the fact that interaction with impurity is local. 
It is integration over k y in Eq. (|B1[) that makes the ex- 
pression to be zero. Indeed, prefactor v y (k) is odd under 
inversion of k y , whereas both Greens functions do not 
change under k y — > —k y . Since this part is present in all 
diagrams in Fig. [3] (b), all these diagrams vanish. 
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